Microscopic calculation of symmetry projected nuclear level densities 
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We present a quantum Monte Carlo method with exact projection on parity and angular momentum that is free 
of sign-problems for seniority-conserving nuclear interactions. This technique allows a microscopic calculation 
of angular momentum and parity projected nuclear level densities. We present results for the 55 Fe, 56 Fe and 
57 Fe isotopes. Signatures of the pairing phase transition are observed in the angular momentum distribution of 
the nuclear level density. 
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The level density is a fundamental property of a many -body 
system as all thermodynamical quantities can be derived from 
it. In nuclear physics, level densities are important because, 
according to Fermi's golden rule, they are critical for estimat- 
ing nuclear reaction rates. The simplest approach to estimat- 
ing the nuclear level density is to assume a free Fermi gas, 
leading to the well-known Bethe formula QjJ, |2fl. A simple 
phenomenological way of incorporating pairing correlations 
and shell effects is to back-shift the excitation energy E x by 
an amount A, resulting in the back-shifted Bethe formula. 

The calculation of statistical nuclear reaction rates also re- 
quires knowledge of the angular momentum distribution of 
the nuclear level density. An empirical formula for the angu- 
lar momentum distribution of the level density at a fixed exci- 
tation energy E x assumes uncorrected and randomly coupled 
single-particle spins, and is given by |H,|3l 



pj(E x ) = p(E x ) 



(2/+1) 
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Here J is the angular momentum and a is a spin-cutoff param- 
eter, which can be used as an energy-dependent fitting param- 
eter. The total level density p (E x ) in (HJ is often parametrized 
by a backshifted Bethe formula. The main drawback of this 
approach is that it requires a fit for each individual nucleus. 

The microscopic calculation of the level density, and in par- 
ticular its angular momentum distribution beyond the Fermi 
gas model and the spin-cutoff parametrization of Eq.([T]i poses 
a complex problem. Recently, progress has been made both 
experimentally [5] and theoretically J6|] in gaining more in- 
sight into the angular momentum distribution of level densi- 
ties. The level density has been calculated microscopically 
using the shell model Monte Carlo (SMMC) method [7, 8]. 
Here, the effect of residual interactions between the nucle- 
ons is taken into account, although restricted to pairing plus 
multipole-multipole interactions that are free of the Monte 
Carlo sign problem. Angular momentum projection in the 
SMMC method was recently carried out for scalar observ- 
ables by first projecting on the angular momentum component 
J z Jh]. Angular momentum projection introduces a new sign 
problem for 7^0 even when for a good-sign interaction. At 



intermediate and high temperatures, this sign problem is suf- 
ficiently moderate to allow the reliable derivation of angular 
momentum distributions. However, at low temperatures this 
sign problem becomes more severe. 

In this Letter, we present a new breakthrough that avoids 
sign problems at any value of angular momentum and/or par- 
ticle number and that allows a microscopic test for the va- 
lidity of the Fermi gas and spin-cutoff model for nuclei that 
exhibit strong pairing correlations. We discuss a quantum 
Monte Carlo (QMC) method to solve a general isovector/ = 
pairing model and solve the angular momentum projection 
problem for seniority conserving models. The pairing model 
serves as a benchmark to identify signatures of pairing cor- 
relations in the nuclear level density. In a finite nucleus, the 
number of paired particles is relatively small, and the identifi- 
cation of such signatures of the nuclear superfluid transition is 
difficult. Including higher order interaction multipoles, would 
also introduce a sign problem in our approach. Nonetheless, 
the projection is direct, and our QMC method allows us to 
treat the very large model spaces, which are needed for the cal- 
culation of level densities at higher excitation energies, and, 
particularly, for the study of the parity and angular momentum 
dependence of level densities (even at low excitation energy). 

We start from a general isovector / = pairing model, here 
written in the quasi-spin formalism ll9i [loll 
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where the operator aj m creates a particle in the spherical 
mean-field single-particle state \ jm) with energy £,, and \ jm) 
is the time-reverse conjugate state of \jm). The operators 
S^SjjSj close an SU(2) algebra and are known as quasi- 
spin operators. The Hamiltonian (|2]i is exactly solvable for 
a constant pairing strength. The QMC method discussed here 
allows us to solve the general case when all terms in <(2J are 
attractive, i.e., gjj > 0. 
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We formulate the pairing problem in the quasi-spin formal- 
ism to address the problem of the angular momentum projec- 
tion. The quasi-spin projection S® determines the total num- 
ber of particles Nj = 25* +£2 ; - in level j = j + 1/2 is 
the pair degeneracy of orbital j), while the quasi-spin quan- 
tum number Sj determines the seniority quantum number 
Vj = £lj — 2Sj, i.e., the number of unpaired particles in level j. 
Within this quasi-spin (or seniority) scheme, angular momen- 
tum remains a good quantum number. Given a set of quasi- 
spin quantum numbers Sj, one can determine the degeneracy 
for a given value of the total angular momentum /. Formally, 
the problem is equivalent to a chain of spins. The pairing inter- 
action can flip spins such that the total "magnetization" Y,j S® 
remains constant (particle number conservation). In addition, 
each quasi-spin quantum number Sj can take values between 
Oand Q.j/2. 

There exist a number of very efficient QMC approaches 
to simulate spin chains I! ill . However, in our quasi-spin 
model, the value of S® is fixed, while the quasi-spin val- 
ues Sj can change. Recently, we developed a non-local loop 
update QMC scheme that is capable of sampling spin mod- 
els at constant magnetization. In our case this means sam- 
pling the canonical ensemble Assuming the Hamiltonian 
H = Ho — V consists of two non-commuting parts Ho and V, 
our scheme is based on a perturbative expansion of the parti- 
tion function at inverse temperature P 

Tr^-P") = £ / dx m / dx*-! ■ 
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with V (t) =exp(— zHo)Vexp(zHo). Here we choose 

Ho = - ^gjjSpj ; V = £ g jf Spj, . (5) 

J J i¥=f 

The basic idea of the QMC method is to insert a so-called 
worm operator A in the partition function, obtaining an ex- 
tended partition function Tr(Ae~$ ). By propagating this 
worm operator through imaginary time according to the rules 
explained in Ref. 11211 . one generates configurations that are 
distributed according to their weight in the canonical parti- 
tion function Tr#(e~' 3 ) at fixed particle number N. Here, 
the worm operator is chosen to be C + Y,ijS^Sj , with C a 
constant. Such a worm operator allows for the sampling of 
all configurations that correspond to a fixed set of quasi-spin 
quantum numbers Sj, without changing the value of Y,j S*j- 
However, to be ergodic, the worm operator must also generate 
configurations with varying seniority quantum numbers, and 
therefore we add a worm operator that can change the values 
of the quasi-spins Si, Sj and their projections 5?, S® for two 
levels i and j such that (Sf + 5^) remains the same (particle- 
number conservation). The complete quasi-spin phase space 
is now sampled by propagating a seniority conserving worm 
operator and an additional seniority non-conserving worm op- 
erator. Additional details on the worm operator moves can be 
found in Ref. d. 



The moves of the worm operator can be constrained in such 
a way that the angular momentum / of the configuration does 
not change. Choosing the initial configuration with well de- 
fined angular momentum quantum number J enables us to cal- 
culate thermal averages for that fixed value of J. Parity is also 
a good quantum number in the quasi-spin basis, so that parity 
projection can be accomplished in a similar way. 

We have used the QMC method outlined above to study 
the angular momentum and parity distribution of nuclear level 
densities in the presence of pairing correlations. We focused 
on the iron isotopes 55 Fe, 56 Fe and 57 Fe, within the complete 
0/ + lp + Qg<) n model space. To study truncation effects, we 
also consider an extended model space Is + Od + Of + lp + 
0g9/2 +2s + Id, The inverse temperature p ranges from to 
2.5 MeV -1 . For the mean-field potential, we used a Woods- 
Saxon potential with the parametrization of Ref. [13]. The 
single particle energies e ; - were obtained by diagonalizing this 
potential in a harmonic oscillator basis. For simplicity, we 
used a constant pairing strength (although our QMC method is 
not restricted to this case), determined to reproduce the experi- 
mental gap parameter - 1/2(® (N- 1 ,Z) - 23 (N,Z) + « (N+ 
1 , Z) ) for Z = 26 protons and N = 30 neutrons, with $ denot- 
ing the binding energy yfl. For the larger valence space, the 
value of pairing strength g is renormalized in such a way that 
the gap parameter remained fixed. Using the temperatures for 
which the heat capacity is maximal, we estimate our calcu- 
lations to be free of truncation effects for excitation energies 
up to E x < 20 Me V and E x < 50 MeV for the p f + g 9/2 and 
sd + pf + g<)/2 + sd valence spaces, respectively. Total and 
angular momentum/parity-projected level densities are calcu- 
lated as usual in the saddle point approximation. 

Figure Q] shows the angular momentum distribution of the 
projected level density pj at four different excitation energies 
E x (we do not include not the magnetic degeneracy in p j, i.e. 
P(E X ) = Ly(2/+ l)pj(E x )). The solid squares are the QMC 
results (the statistical error are much smaller than the size of 
the squares), while the solid lines are the spin-cutoff model 
Eq. (Q]) fitted to the QMC data with a 2 as a fit parameter. For 
all three iron isotopes, the angular momentum distribution be- 
comes broader with increasing excitation energy. When this 
energy is sufficiently high (E x > 10 MeV), the QMC data are 
well described by the spin-cutoff model. For excitation en- 
ergies < 10 MeV, some deviations are observed. In particu- 
lar, a staggering effect (in J) is found for the even-even 56 Fe 
isotope but not in the odd-even Fe isotopes. This effect was 
also reported in Ref. (§], where the interaction included higher 
multipoles, and was recognized as a signature of pairing cor- 
relations. Here, we observe a strong suppression of the J = I 
level density (pi « po), which is much smaller than the J = 2 
level density. For higher angular momenta, we found only a 
very small staggering in the angular momentum dependence. 
This can be understood from the fact that the pairing interac- 
tion scatters only 7 = nucleon pairs. 

In a semiclassical approach, the spin-cutoff parameter is 
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FIG. 1 : The angular momentum distribution of the level density at a 
given excitation energy E x for the iron isotopes 55 Fe, 56 Fe and 57 Fe. 
Solid squares arethe Monte Carlo results in the pf + ggn valence 
space. The solid lines are fits to the spin-cutoff model Eq. 10. Also 
shown are distributions given by the spin-cutoff model with c 2 cal- 
culated from Eq. (f7) using the rigid body moment of inertia (dotted- 
dashed lines), and from Eq. l|6j (dotted lines). 



given by 

a 2 = ^>, (6) 

with (. . .) the thermal average at temperature T lfl4ll . It is then 
common to define an effective moment of inertia from 



The spin-cutoff parameter is often determined by substitut- 
ing the rigid-body moment of inertia, / = 2mA(roA 1 / 3 ) 2 /5 in 
Eq. Ql. (Here, ro is the nuclear radius parameter, A is the 
mass number and m is the nucleon mass). The spin-cutoff 
model (0 with o 2 calculated from the rigid body moment of 
inertia (dotted-dashed lines in Fig. [T} essentially coincides for 
E x > 10 MeV with the distribution fitted to the Monte Carlo 
data (solid lines in Fig. [TJ. For 56 Fe, however, for E x < 10 
MeV, the rigid-body moment of inertia predicts a broader an- 
gular momentum distribution than the one described by the fit. 
This indicates a reduction of the effective moment of inertia 
because of pairing correlations, and is a signature of nuclear 
superfluidity. We also determined o 2 from the thermal average 
(J 2 ), calculated directly in the QMC simulation. The angular 
momentum distributions using these spin-cutoff parameter are 
shown in Fig. [TJby the dotted lines. These distributions are 
slightly more peaked than the fitted distributions (solid lines), 
especially at low excitation energies (E x < 10MeV).For 56 Fe 
at E x = 5.5 MeV, this curve coincides with the fit. 



The top panels of Fig. [2]shows the effective moment of in- 
ertia as a function of excitation energy. At high excitation en- 
ergy, the moment of inertia calculated from (J 2 ) is very close 
to its rigid body value, the latter indicated with the horizontal 
dashed lines. We also calculated the spin-cutoff parameter o 2 
directly from the ratio po(i/2)/P m the even (odd) case. The 
corresponding moments of inertia are also shown in the top 
panels of Fig. [2] (solid circles). These values seem to differ 
from our previous estimates of the moment of inertia, but they 
do not lead to significant differences in the angular momen- 
tum distributions, except for 56 Fe for E x < 10 MeV. Here, the 
moment of inertia is strongly suppressed compared to 55 Fe 
and 57 Fe, due to J = pairing. 

The bottom panels of Fig. [2] show the pair correlation en- 
ergy (S + S~) = (Ej.j' SjSji) as a function of excitation en- 
ergy. It is seen that the pair correlation energy is strongly 
suppressed with increasing excitation energy. The suppressed 
neutron correlation energy of 55 Fe and 57 Fe is seen to increase 
first slightly with increased excitation energy (or temperature). 
This results from a blocking effect of the unpaired neutron at 
low temperatures. 
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FIG. 2: Top panels: the moment of inertia / (in units ft 2 ) for 55 Fe 
(left), 56 Fe (middle) and 57 Fe (right). The horizontal dashed lines 
are the rigid body moments of inertia, the solid lines are moments of 
inertia calculated from Eq. and the solid circles are from Eq. ((7} 
with o 2 determined via Po/p (in the even-even case) and P1/2/P ( m 
the odd-even case). Bottom panels: the pair correlation energy (S'S). 



In Fig. [3] we show the parity-projected level density for 
/ = 0,1,2,3 as a function of excitation energy using the 
sd + fpg + sd valence space. For energies below ~ 20 MeV, 
the angular momentum projected level densities show a strong 
parity dependence. Recently, angular momentum and parity- 
projected level densities were determined experimentally |5J]. 
However, no parity dependence was resolved for level densi- 
ties with specific values of J. Our calculations show that the 
angular momentum projected level density displays a strong 
parity dependence. For J = 2, we also show the level density 
within the smaller pf + gg/2 valence space (dashed lines). Be- 
low ~ 20 MeV, the even-parity level density is in good agree- 
ment with the results found in the sd + fpg + sd model space. 
However, the inclusion of the sd shells, below and above the 
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Pf + 89/2 shell, significantly enhances the odd-parity level 
density at low excitations. We found a similar effect for the 
total odd-parity level density. This results from the increased 
fraction of single-particle levels with positive parity. 
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FIG. 3: The angular momentum and parity-projected level density of 
56 Fe for the four lowest angular momentum values. A strong parity 
dependence is found for excitation energies up to ~ 20 MeV. The 
valence model space is sd + pf + gg/2 + sd- For J = 2, we also show 
results for the smaller pf + gg/2 valence space (dashed lines). 

In the left (right) panel of Fig. @] we show the QMC to- 
tal level density, as well as the angular momentum projected 
level densities 7 = (7 = 1) and 7 = 2 (7 = 3) for 56 Fe. 
The QMC total level density is well described by the back- 
shifted Bethe formula (fitted for 4.5 MeV < E x < 40 MeV) 
with Fermi gas parameter a — 5.741 ±0.034 MeV -1 and 
backshift A = 1.591 ±0.057 MeV (see dotted-dashed line in 
Fig- ID- These are similar to the SMMC level density parame- 
ters a = 5. 780 ±0.055 MeV" 1 and A = 1.560±0.161 MeV of 
Ref. (Sfl, which were found to be in good agreement with the 
experimental level density. The present results indicate that 
the backshifted Bethe formula with temperature-independent 
parameter a works well up to higher energies [15]. The dashed 
lines of Fig. 0]are the angular momentum projected level den- 
sities, calculated from the total level density through Eq. (Q} 
with rigid body o 2 . At high excitation energy, these densities 
coincide with the QMC densities. However, for excitation en- 
ergies below ~ 10 MeV, the QMC 7 = level density deviates 
from the level density predicted by the spin-cutoff model with 
rigid body a 2 . This is a signature of the pairing phase transi- 
tion in the level density. This signature is also visible in the 
7 = 2 level density for energies below ~ 8 MeV. For the odd 7 
values, however, all QMC data are well described by the rigid 
body model. 

In conclusion, we have used a quantum Monte Carlo 
method to calculate angular momentum and parity-projected 
level densities for a nuclear pairing model. This method al- 
lows us to calculate level densities in the very large model 
spaces that are necessary to avoid truncation effects. We have 
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FIG. 4: The QMC total level density for 56 Fe, well described by 
the backshifted Bethe formula (dotted-dashed line), together with the 
projected densities J = 0, J = 2 (left panel), and J = 1,7 = 3 (right 
panel). The solid lines connect the 7-projected Monte Carlo results, 
while the dashed lines show the projected level density from the spin- 
cutoff model Eq. |QJ using the rigid body values of o 2 . 



found that pairing correlations affect the angular momentum 
distribution of the level density at low excitation energies, 
thereby revealing signatures of the pairing phase transition in 
nuclei. 
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